Mon. Not. R. Astron. Soc. 000, 000-000 (1994) Printed 1 February 2008 (MN plain macros vl.6) 



How to drive roAp stars 



Alfred Gautschy^, Hideyuki Saio^, Housi Harzenmoser^, 

i Astronomisches Institut der Universitdt Basel, Venusstr.7, CH-4102 Binningen, Switzerland 
Astronomical Institute, Tohoku University, Sendai, Japan 

ABSTRACT 

Rapidly oscillating Ap stars constitute a unique class of pulsators to study nonradial 
oscillations under some - even for stars - unusual physical conditions. These stars are 
chemically peculiar, they have strong magnetic fields, and they often pulsate in several 
high-order acoustic modes simultaneously. We discuss here an excitation mechanism 
for short-period oscillation modes based on the classical k mechanism. We particularly 
stress the conditions that must be fulfilled for successful driving. Specifically, we discuss 
the roles of the chemical peculiarity and strong magnetic field on the oscillation modes 
and what separates these pulsators from S Scuti and Am- type stars. 
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1 CHARACTERIZING THE BREED 

Rapidly oscillating Ap (roAp) stars constitute a low- 
temperature subgroup of the chemically peculiar A-type 
(Ap) stars with strong Sr, Cr, and Eu overabundances. The 
, first star of this class, HD 101065, was discovered photomet- 

■ rically in 1978 by Kurtz (Kurtz 1978, Kurtz & Wegner 1979). 
Presently, the group counts about 30 members. The periods 

I of the light variability, which can be multi-periodic, range 
from roughly 4 to 16 minutes. As the roAp stars lie in the 
instability strip close to the main-sequence they are thought 

. to have masses between 1.5 and 2Mq. Therefore, the short 

' periods must be interpreted as high radial-order p modes. 

, The amplitudes of variability in the blue, where they are 

■ largest, stay below about a hundredth of a magnitude. 

The roAp stars show strong surface magnetic fields 
, exceeding frequently some kilo Gauss (cf. Mathys & Hu- 
' big 1997) with the magnetic axis being inclined relative to 
the star's rotation axis. The rotation with periods above 

2 days leads, therefore, to a modulation of the magnetic 
field strength. This property gave rise to the oblique rota- 
tor model (Stibbs 1950). From the temporal phase behavior 
of the magnetic field during a rotation period, a dominant 
dipole topology was derived for these stars' global magnetic 
fields. 

The observed chemical peculiarities are interpreted as 
being a very superficial phenomenon which is restricted to 
the upper stellar atmosphere. It is commonly believed that 
the strong magnetic field suppresses turbulent motion which 
otherwise quickly destroys any diffusively built up segrega- 
tion of atomic species that react differently on the prevailing 
stellar radiation field (e.g. Michaud 1980, Vauclair & Vau- 
clair 1982, Alecian 1986) . None of the available diffusion the- 
ories reproduces, however, quantitatively the observed vari- 
ety of peculiar-star characteristics with a satisfactory low 
number of freely tunable parameters. 

The short-period oscillation amplitudes of the roAp 
stars are temporally not stable but they vary cyclically in 



phase with the magnetic field modulation. The variability 
is in accordance with the postulation of dipole oscillation- 
modes being aligned with the magnetic axis of the roAp stars 
which changes its aspect angle during the star's rotation pe- 
riod. This picture - called the oblique pulsator model - is 
furthermore supported by phase jumps of the pulsational 
phase at quadrature of the magnetic axis with respect to 
the line-of-sight. If the the temporal variation of the pulsa- 
tion is parameterized by cos(ijJt-\-ip), then the phase ifi jumps 
by TT radians at the moment of magnetic polarity reversal; 
and this is just what is observed (cf. Kurtz 1990) . 

The roAp variability is unique in several respects. De- 
spite these stars sharing their territory in the instability strip 
with 5 Set stars, they oscillate with much shorter periods. 
Much higher overtones are obviously excited in roAp stars. 
The mode density in the corresponding frequency domain 
is high, and it remains unclear what selection mechanism 
is picking out the very low number of observed oscillation 
modes. The spherical degrees of the modes identified in some 
cases (e.g. Kurtz 1990) were axisymmetric dipole modes. 
In general, they must definitely all be axisymmetric, low-£ 
modes. 

Rotation and magnetic field of the roAp stars cause 
the m-degeneracy of the oscillation frequencies to be lifted. 
From the magnitude of the frequency splitting and the am- 
plitudes of the multiplet components, aspect angles of the 
rotation and magnetic axis and spatially integrated mag- 
netic field strengths can be deduced. If several radial or- 
ders of oscillation modes are observable (as e.g. in HR 1217, 
Kurtz et al. 1985) frequency spacings can be deduced. As 
roAp stars oscillate in high radial orders and low spheri- 
cal degree, asymptotic pulsation theory, together with other 
(spectroscopic, photometric, astrometric) observables might 
constrain global stellar parameters such as mass and age (cf. 
Christensen-Dalsgaard 1988, Brown et al. 1994). 

The lifetimes of oscillation modes in roAp stars range 
from days to years. For long-lived modes, cyclically (but not 
periodic) varying oscillation frequencies on the timescale of 
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months were reported recently (Kurtz et al. 1994, Kurtz 
et al. 1997). The origin of these shifts is unclear. Specula- 
tions, however, tend towards the action of magnetic activity- 
cycles as soon on tho sun. All in all, tlio roAp stars con- 
stitute a unique laboratory to test important domains of 
stellar pulsation theory and of stellar-hydromagnetics com- 
municated through oscillations. For more details we refer 
to recent reviews which addressed comprehensively most - 
theoretical as woU as observational - aspects of roAp stars: 
Kurtz (1990), Shibahashi (1991), Matthews (1997). 

In Section 2 we address the stellar evolution aspect of 
the roAp stars and introduce the models with which we 
worked in our analyses. Section 3 deals with the nonadia- 
batic oscillation spectra of the roAp stars and tho possible 
underlying excitation mechanism. Success and limitations of 
our work are critically evaluated and compared with results 
in the literature in Sect. 4. Finally, we devote an Appendix 
to the comparison of adiabatic and nonadiabatic oscillation 
spectra. 



2 ASPECTS OF STRUCTURE AND 
EVOLUTION 

The overlapping of the S Set variables and of roAp stars 
in the instability strip on the HR diagram suggests that the 
latter have similar masses and that they arc also in a compa- 
rable evolutionary stage. This is additionally supported by 
the asymptotic frequency separations (cf. Appendix) which, 
for d Set-like stars, are compatible with those seen in roAp 
stars (e.g. Gabriel et al. 1985, Shibahashi & Saio 1985, Heller 
& Kawalor 1988). Therefore, wc analyzed stellar models - 
representing roAp stars - with masses between 1.5 and 1.87 
Mq during the main-sequence and the early hydrogen shell- 
burning evolutionary phase (cf. Fig. 1). 

Referring to Hipparcos observations. North et al. (1997) 
attempted recently to pin down the location of the roAp 
stars on the HR diagram. The result is plotted in Fig. 1 
with grey shading. The dark grey area enclosed the mean 
values obtained for the sample of 14 stars. The lighter grey 
encloses the max;imum extensions of the error bars specified 
in North et al. (1997). Their data suggest that stars with 
masses above roughly 1.6 Mq could become roAp stars un- 
der suitable conditions mainly during their main-sequence 
stage. The dotted line in Fig. 1 indicates the position of 
the observed blue edge for 5 Set variables. It is evident that 
an overlap exists; therefore, roAp and S Set stars appear to 
share common ground. 

For reference, the locus of vq = const, is plotted on the 
HR diagram in Fig. 1. The mathematical definition of i^o - 
the asymptotic frequency spacing of equal-degree modes - is 
given in the Appendix. The determination of this observable, 
together with a reliable temperature measure, constitutes an 
important characterization of a star. The z/q values of our 
models are listed in Tab. 1. 

The stellar evolution computations for this study were 
performed with the Basel stellar evolution code using the 
micro-physics described in Gautschy et al. (1996). The ho- 
mogeneous chemical composition of the stellar models on the 
ZAMS was in all cases X = 0.7, Y = 0.28. When evolving 
the stars we did not account for any chemical peculiarity. We 
assumed that an atmosphere with a prescribed T—t relation 
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Figure 1. Hcrtzsprung- Russell diagram with the loci of the early 
tracks of main-scqucncc and early post-main sequence evolution 
of stars with 1.5, 1.6, 1.7, and 1.87 Mq. The positions along the 
tracks where stellar models were analyzed for their oscillation 
behavior are marked with capital letters. The corresponding pa- 
rameters are listed in Tab. 1. The grey areas denote the map of 
observed roAp stars onto the HR plane as published by North et 
al. (1997). For details see the text. The dotted line indicates the 
location of the observed blue edge of S Set stars. 
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Figure 2. Logarithm of /3, the ratio of magnetic pressure to total 
pressure throughout a star in the direction of the dipole mag- 
netic pole. The surface magnetic field strength is chosen to be 
10'^ Gauss. The stellar model corresponds to model B which is 
characterized in Tab. 1. 

is overlying the stellar photosphere where the stellar struc- 
ture computations stopped. Any changes in chemical com- 
position and molecular weight which characterize Ap-star 
peculiarity were restricted to these atmosphere calculations. 
Descriptions and results of specific choices arc presented in 
Sect. 3. Furthermore, we did not consider effects of a dipole 
magnetic field in the computation of the stellar structure. 
In this respect it is worthwhile to notice that even a kilo- 
Gauss non-potential field is hardly of relevance for most of 
the star: Figure 2 shows the logarithm of the ratio of the 
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Table 1. Stellcir parameters of models used for stability analyses. 

10'^ X 3GM/R^ 
5.626 
1.611 
1.036 
5.917 
2.257 
1.349 
5.967 
2.308 
1.365 
5.712 
2.217 
1.348 



Model 


log Teff 


log L/Ls 


M/Mq 




A 


3.936 


1.086 


1.87 


88.36 


B 


3.884 


1.240 


1.87 


51.46 


C 


3.880 


1.352 


1.87 


.39.37 


D 


3.904 


0.916 


1.70 


89.89 


E 


3.864 


1.035 


1.70 


56.56 


F 


3.864 


1.170 


1.70 


44.26 


G 


3.880 


0.804 


1.60 


90.30 


H 


3.843 


0.927 


1.60 


57.46 


I 


3.843 


1.079 


1.60 


44.77 


J 


3.853 


0.686 


1.50 


90.79 


K 


3.816 


0.812 


1.50 


57.67 


L 


3.816 


0.960 


1.50 


45.34 



magnetic pressure to the total pressure at the magnetic pole 
(for a dipole field) for stellar model B. Evidently, the mag- 
netic pressure exceeds the hydrostatic pressure in the very 
outermost regions only. Since the the stellar envelopes of 
the stellar models corresponding to roAp stars are mainly 
radiative they approach quickly the radiative-zero solution. 
Therefore, surface information is not communicated to the 
deep interior which drives the star's evolution. The run of 
(5 in Fig. 2 stops at the outer edge of the convective central 
region where the the dipole nature probably breaks down. 

Ftom the stellar evolutionary sequences we selected a 
few models on which the stellar oscillation computations 
were performed. For each of the tracks - 1.5, 1.6, 1.7, and 
1.87 Mq - we selected two models each in the core hydrogen- 
burning phase and one in the early hydrogen shell-burning 
stage as representatives of roAp stars. Table 1 lists, for mod- 
els A to L, some relevant stellar and oscillatory parameters. 



3 THE RAPID OSCILLATIONS OF Ap STARS 

In this section we focus on pulsational aspects of roAp stars. 

First, we address their oscillation-frequency spectra and the 
problems they pose in understanding the underlying stellar 
models. The later part of Sect. 3 is devoted to nonadiabatic 
aspects. In particular we suggest an excitation mechanism 
for the observed high-order p modes. 

The frequently occurring dimensionl ess frequencies a = 
(o"R,cri) are normalized by -y/SGMi/R-J. All the numerical 
computations of oscillation properties were performed with 
codes similar to the one described in Gautschy et al. (1996). 



3.1 Oscillation-mode physics 

To get an impression of the mode cavities in stars it is appro- 
priate to plot the characteristic frequencies ctc- They sepa- 
rate the sjjatially propagative from the evanescent regions 
in a star (see e.g. Unno et al. 1989). We compute the disper- 
sion relation of the adiabatic Cowling equations assuming a 
plane- wave with a spatial behavior of exp(ifc • x) in any per- 
turbed physical quantity, where we use x = r/R*. Searching 
for the loci of Real(/c) = leads to a quadratic for ctc- 
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The quantity A* = — V^(l/ri — d log p/d log P) is essentially 
the Schwarzschild determinant. For the meaning of the re- 
maining symbols we refer to Unno et al. (1989). 

As an example of the propagation properties of dipole 
modes in roAp-like stars we show in Fig. 3 the characteristic 
frequencies of stellar model B which is our slightly evolved 
1.87M0 star. The upper and lower lines represent ac+ and 
ac-, respectively. An oscillation mode with frequency a 
propagates spatially as a p mode in the regions where cr 
exceeds the larger value of ac- The dominant features of 
the upper ac line are the two depressions in the superficial 
regions which are caused by partial ionization of H and He 
(4 < logP < 4.7) and He+ (5.5 < logP < 5.8). The top 
of the H/He ionization zone reaches into the photosphere. 
Furthermore, the steep composition gradient which builds 
up at the outer edge of the nuclear burning core causes the 
glitch in ac at about log P = 16.8. The influence on the 
cavities of the composition change is even more pronounced 
in the lower ac curve which traces the upper boundary for 
the propagation region of the g modes. 

The oscillation frequency domain observed in roAp stars 
is marked by the grey area shown in Fig. 3. We see that 
all the modes in the appropriate frequency range are pure 
acoustic modes. A slight influence on the mode spectrum 
from the steep molecular-weight (/i) gradient at the outer 
edge of the nuclear burning convective core can be expected. 

After the roAp variability was interpreted to be due 
to high-order acoustic oscillation modes, the magnitude of 
the critical frequency (ac) in the stellar atmosphere, i.e. at 
the outer boundary of the oscillation cavity, became a point 
of concern. The superficial acoustic cut-off frequency, aac, 
which corresponds to ac there, determines if a wave with 
frequency a is refiected (if an < fac) or if it is partially trans- 
mitted (if CTac < ctr) and therefore loses part of its energy 
into the atmosphere. Stellar models which are appropriate 
for roAp stars show the convection zone induced by H/He 
ionization to reach up to the stellar photosphere. At the pho- 
tosphere the acoustic wave propagation is determined by the 
run of the Brunt- Vaisala frequency which is imaginary in the 
convectively unstable layers. To find the outer edge of the 
acoustic cavity, some atmospheric modeling is necessary. In 
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Figure 3. Characteristic frequencies of dipole modes for model 
B as function of depth which is measured in log P. The grey-filled 
area covers the observed oscillation-frequency domain of roAp 
stars. 

other words, the action of the outer boundary condition is 
dominated by the pecuMarities in the optically thin radia- 
tive layers. The use of standard (e.g. Eddington-grey) T — t 
relations leads to rather low cut-off frequencies. Many of the 
observed oscillation frequencies lie above the cut-off so that 
some mode energy is lost from the star (cf. Fig. 3). If the 
cut-off frequency in roAp stars were indeed that low, then 
the excitation mechanism would have to be very efficient to 
over-compensate the modal leaJs into the atmosphere. 

To resolve the dilemma of the low cut-off frequen- 
cies at the surfaces, considerable efforts went into hypoth- 
esizing physical processes which act in the stellar atmo- 
spheres to sufficiently increase the acoustic cut-off frequency. 
Shibahashi & Saio (1985) computed critical frequencies for 
Eddington-type atmospheres as well as for T — t relations 
fitted to more detailed Kurucz atmospheres. Their fitted so- 
lutions showed a larger ratio of Tcs/T{t — 0) than in the 
Eddington atmospheres. For none of the chosen T — t func- 
tions, however, a fXac resulted that would lead to fully re- 
flected roAp-relevant oscillation modes. The results of an 
observational determination of the T — t relation of a roAp 
star were presented by Matthews et al. (1996). Their studied 
star, HR 3831, indicated am even steeper temperature drop 
close to the photosphere and a more isothermal behavior 
at small optical depths than what Shibahashi & Saio (1985) 
suggested. (It is, however, noteworthy that Kurtz & Mcdupo 
(1996) in their analysis of frequency-dependent oscillation 
amplitudes conclude that a T — r relationship cannot be de- 
rived from roAp oscillation data.) A list of theoretical sug- 
gestions to increase (Jac at the outer boundary of the stellar 
models was given in Section 5 of Shibahashi (1991). 

To simulate the effect of roAp-type chemical inhomo- 
geneities on the characteristic frequencies, in particular on 
the atmospheric cut-off frequency, we studied sharp /i transi- 
tions in atmospheres with standard T — t relations. Accord- 
ing to eq. 3, the floating on the top of the stellar atmosphere 
of rare earth elements leads to a fj,- inversion. The magnitude 
and the location of such //-jumps were guided by results of 
Michaud (1980). Our atmospheric /i variations led to ac de- 
pressions which had no noticeable influence on the amount 
of energy lealcage of the waves with frequencies above the 
cut-off frequency. Acoustically, the feature was found to be 
too shallow to be noticed even by high-frequency p modes. 

We postulate now a mechanism which rather success- 
fully increased ac in our models. It remains to be seen, 
however, if it is indeed realized in nature: Since roAp stars 
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Figure 4. Dependence of the characteristic frequency ctq (bot- 
tom panel) in the superficial layers on the choice of the T — t 
relation (top panel) in the star's atmosphere. The case of dipole 
modes in model B are depicted. 

have a convection zone which reaches the photosphere and 
since they have strong magnetic flelds, we suggest that roAp 
stars are capable of building up - very much like solar-type 
stars - a chromosphere of some extent. This chromosphere 
goes along with a temperature inversion which we ad hoc 
parameterized by: 

T(r) = 0.931 • Tcfi • [r + 0.72858 - 0.3367 • exp(-2.54r) 

+ 10 • exp(-50r) - 1.0 ■ exp(-500r)]^''''. 

The above choice is motivated by the Shibahashi & Saio 
(1985) fit (full line in Fig. 4) except that the first coefficient 
in the square bracket is adapted to enforce T(r — 2/3) = 
Toff . The positive sign in front of the exp(— 50r) coefficient 
of cq. 2 induces a temperature inversion at small optical 
depths (dash-dotted line in Fig. 4). Replacing this term by 
-|-120-exp(— 40t) leads to the stronger temperature inversion 
delineated with the dashed line. Finally, the dotted line in 
Fig. 4 indicates the run of an Eddington-type T — t relation. 

We cam understand the effect of the temperature inver- 
sion on the p-mode cavity by approximating ac for acoustic 
modes close to the surface: 



ac 



V ( dlogp 1 



3 V^dlogP Fi 
V 



(2) 



(3) 



We used a = (Slog p/9 log P)t, & = —{d\og p/ dXogT)^, 
and Vt = dlogT/dlogP to be the effective temperature 
gradient. A negative Vt clearly increases the critical fre- 
quency for adiabatic p modes. Quantitative effects are dis- 
played in Fig. 4. 

The temperature inversions shown in Fig. 4 are, admit- 
tedly, rather strong. In particular, we have to face the fact 
that no observational signature for chromospheric activity 
has been found to date (e.g. Shore et al. 1987). Therefore, 
the conjectured temperature inversion cannot be very strong 
because it should have been observable through emission line 
features otherwise. Notice, however, that only the innermost 
part of the above temperature inversions are essential for the 
increase of ac - Hence, a leveling off of the temperature after 
an increase of 3 000 Kelvin at most is sufficient to completely 
reflect the roAp-type acoustic waves in many cases. 
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In the next section we investigate the effect of this mod- 
ified outer edge of the cavity on the stability properties of 
high-order p modes. As of now, the physics of their excita- 
tion remains a matter of debate. 

3.2 Excitation physics 

Since the roAp star domain on the HR diagram appears to 
overlap with the 5 Set stars (see Fig. 1), the partial He"*" 
ionization zone was thought for a long time to be also the 
driving agent for roAp stars (e.g. Kurtz 1990, Matthews 
1997). Numerical computations never supported this con- 
jecture, however. Based on the compositional peculiarity 
of the roAp stars, Matthews (1988) hypothesized, relying 
on qualitative micro-physical arguments, that Si IV par- 
tial ionization might provide sufficient p-mode driving even 
if He was drained out of the potential driving depth in 
the star. He estimated an overabundance of Si IV in the 
superficial stellar regions of the order of 200 to be suffi- 
cient for the mechanism to work. A completely different av- 
enue was taken by Shibahashi (1983) which was also inves- 
tigated by Cox (1984). Their local analyses of oscillatory 
magneto-gravity waves suggested overstability under favor- 
able circumstances. They expected the magnetic polar re- 
gions to become overstable under Ap-star circumstances. 
The eventually observable roAp star oscillation are those 
global p modes which are resonantly excited by the over- 
stable magneto-gravity waves. This picture has the great 
advantage of naturally explaining the alignment of the pul- 
sation axis with the magnetic axis and the selection of ax- 
isymmetric modes. Detailed global analyses which are nec- 
essary to settle the excitation problem eventually are, how- 
ever, lacking also for this approach. Recently, Dziembowski 
& Goode (1996) noted that the partial H ionization con- 
tributes positively to the work integral. No p mode with a 
roAp-star relevant period was found yet destabilized by the 
K mechanism. 

If a "thermo-mechanical valve mechanism" is to power 
roAp-like oscillation modes, then the oscillation periods have 
to be of the order of the thermal timescale (rtherm) of the 
envelope overlying the driving region (e.g. Cox 1980). The 
corresponding spatial structure shown in Fig. 5 for stel- 
lar model E is representative for the whole model series 
looked at in this paper. A qualitative inspection of the 
models shows that those superficial regions with a thermal 
timescale below about 1 000 s contain the steep kt rise to- 
wards the surface induced by partial H and He ionization 
(see Fig. 5). The deeper- lying, potentially driving region for 
which dKx/dr > (at 4.1 ;g logT ^ 4.8 in Fig. 5) favors 
time-scales - as expected - of modes in the 5 Set period 
domain, i.e. above 10'' s. Therefore, we expect - already a 
priori - partial ionization of H and He to be the prospective 
excitation agents for roAp-type oscillations. 

We computed nonadiabatic oscillation spectra for the 
stellar models listed in Tab. 1. The chromosphere-like tem- 
perature structure in the atmosphere which we introduced 
provided a sufficiently high acoustic cut-off so that the 
acoustic cavity appears refiective at the surface for most 
of the computed oscillation modes. The numbers of over- 
stable p modes in the different models are listed in Tab. 2. 
For the 5 Set-type oscillation modes we indicate only if we 
found any ("x" in Tab. 2) or none ("-" in Tab. 2). Selected, 
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Figure 5. Spatial run of kt, a scaled k, and (in the lower panel) 
the variation of the thermal time-scale (rtiierm) in model J. The 
grey shading indicates the temporal region of interest for roAp 
oscillations (cf. text). 
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Figure 6. Appropriately scaled imaginary parts of p modes for 
the 3 models of the I.6M0 sequence with the weak tempera- 
ture inversion. Full line: model G, dotted line: model H, and 
dashed-dotted line: model I. 

results are discussed in the following to demonstrate their 
properties and the physical mechanisms. 

Figure 6 shows the whole dipole p-mode domain in 
which overstable oscillation modes were found for the I.6M0 
models. The imaginary parts of the eigenfrequencies were 
suitably scaled to accommodate the large numerical range 
of cTi in the short- and long-period domain. We defined: 



Ei = Sign (o-i) • log 1 + 



1 ■ 10- 



(4) 



Negative Ei values represent overstable modes. Values of 
Ei exceeding unity behave essentially logarithmic whereas 
smaller values are linear in ai. 

The full line in Fig. 6 (as well as in other figures show- 
ing the three evolutionary stages of a selected stellar mass) 
stands for the zero-age main-sequence (ZAMS) model, the 
dotted line represents the slightly evolved model close to the 
first turn-around, and the dashed-dotted line depicts the re- 
sults of the most evolved stage. For all models in Fig. 6 
we used the weak temperature inversion. The I.&Mq ZAMS 
model is pulsationally stable over the whole frequency do- 
main (i.e. roAp and 5 Set modes). The more evolved models 
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Table 2. Excited dipole modes of roAp and S Set type for weak (first entry) and strong (second entry) T- inversions and 
homogeneous envelope compositions. 



Period domain 


A 


B 


C 


D 


E 


F 


G 


H 


I 


J 


K 


L 


roAp 


0,0 


0,0 


0,0 


0,0 


0,0 


0,0 


0,0 


3,3 


3,2 


1,3 


5,4 


4,3 


5 Set 


~5~ 


x,x 


xr 


x,x 


x,x 




X,- 


X,- 


X,- 


x,x 


x,x 


x,x 



show a short-period instability around 500 s (the number of 
ovorstablo modes is listed in Tab. 2) as well as overstable 
long-period modes above 3 600 s, which are appropriate for 
6 Set variables. For both stellar models (H and I) the range 
of overstable roAp-type modes remains unchanged. This is 
different in the case of a strong T-inversion in the 1.5Mq se- 
quence which is shown in Fig. 7. In this case the instability 
domain shifts to longer periods as the star evolves. On the 
ZAMS, periods around 400 s are excited where the subgiant 
model has 4 overstable p modes around 750 s. 

In all models with overstable high-order p modes we 
found only a few of them (cf. Tab. 2) at a selected spherical 
degree. This short period instability domain is always sepa- 
rated from the long-period domain with low-order p modes 
which are typical for 5 Set variability. It is only in the 1.5Mq 
ZAMS model where even this longer-period instability range 
extends to below 1 000 s. 

To study the destabilisation of the p modes we plotted 
some relevant physical quantities of model H as a function of 
depth in the envelope (parameterized by the total pressure) 
in Fig. 8. The top panel shows the opacity derivatives along 
the chosen thermodynamic basis. It was argued in Unno et 
al. (1989) that a necessary condition for pulsational driving 
requires 



_d_ 

dr 



KT + 



Fa 



> 0. 



(5) 



The major contributor to the above sum is clearly kt- The 
outer potential driving region is due to partial H and He ion- 
ization. The second, deeper lying one is associated with the 
He^ partial ionization. We see that the long-period mode is 
mainly driven by He'*' ionization, the outer ionization zone 
is actually damping in this case. The high-order {k = 34) 
dipole mode, on the other hand gets all the driving from 
the H/Hc ionization zone. The deeper lying regions are all 
slightly damping. The results displayed in Fig. 8 agree very 
well with the rather qualitative reasoning on the potential 
driving agent for roAp modes used before. Notice, how- 
ever, that the condition in Eq. 5 was derived for weakly 
non-adiabatic circumstances. Therefore, noticeable devia- 
tions from this simple form must be expected occasionally. 

3.3 Inhomogeneous envelopes 

A persistent problem with the classical k instability mech- 
anism found above is that not only short-period modes of 
high radial order are excited, but also the 6 Set-type low- 
orders. Observations clearly show that S Sct-Iike oscillations 
are not present in roAp stars. 

To suppress the 5 Set-type modes we assumed a 
drainage of He from the superficial regions. Stabilisation of 
turbulent bulk motion by the magnetic field might cause 
such a He depletion in the outer parts of Ap-star envelopes 
(Vauclair & Vauclair 1982). We assumed ad hoc that the He 
drainage extends over those layers where magnetic pressure 
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Figure 7. Appropriately sealed imaginary parts of p modes 
for the 3 models of the 1.5Mq sequence and a strong temper- 
ature inversion. Full line: model J, dotted line: model K, and 

dashed-dotted line: model L. 
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Figure 8. Top panel shows the spatial run of the opacity deriva- 
tives which enter the computation of the work integral. The 
lower panel displays the cumulative work done by a short-period 
(dashed line) and a long-period (full line) mode, both of which 
are overstable. 

(for a kilo-Gauss field) exceeds gas pressure. We depleted he- 
lium to a mass fraction of 0.1 (and enhanced H to 0.88) at 
temperatures below logT = 4.5 (cf. Fig. 2). Across a tran- 
sition region of AlogT = 0.1 our canonical population I 
composition of X = 0.7, Y = 0.28 was recovered. For ease of 
computation, the smooth transition was prescribed by a hy- 
perbolic tangent function. Fig. 2 indicated already that only 
an extremely thin mass layer is affected by the drainage. In 
models H, K, and L for which we performed this experi- 
ment the He drainage extended over the outer 4 • 10~®M, . 
Notice that this is six orders of magnitude less than what 
Shibahashi & Saio (1985) adopted in their study. 

The repeated nonadiabatic stability calculations on 
models H, K, and L showed that the long-period instabili- 
ties disappear if the superficial layers are chemically inho- 
mogeneous. The effect is less pronounced in model H than 
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Figure 9. The run of Si as a function of period for model K 
with strong temperature inversion. The full line shows the results 
for the inhomogeneous envelope. The dashed line depicts the be- 
havior of oscillation modes in the model with homogeneous outer 
layers. 

in tiie two 1.5Mq ones, iiowcvor. Additionally, tlic oscilla- 
tion modes of tiie strong temperature- inversion models seem 
to react more strongly to inliomogeneity than tiiosc witli 
the weak inversions. The dipole-mode results of model K 
are shown in Fig. 9. The long-period modes above about 
1 200 s which are overstal)le in the model witli cliemically 
homogeneous envelope (broken line) are all stabilized in the 
He-drained one (full line). The short-period modes remain 
overstable. The radial orders of the overstable modes are, 
however, lower in the inhomogeneous envelope: In the ho- 
mogeneous outer layers, periods around 500 s are excited 
whereas periods of about 750 s are favored in the inhomoge- 
neous case. This property is not generic, however. In other 
cases, such as model L, the short-period instability domain 
remained unchanged. A peculiarity of this model L was that 
in addition to the short-period iustability a dipole-mode at 
about 900 s turned overstable in the inhomogeneous enve- 
lope model. In other words, we had one case where non- 
successive radial orders were excited. We show the corre- 
sponding diagram in Fig. 10. 

To see where the loss of the 5 Set modes in inhomo- 
geneous envelope models comes from, we look at the work 
integral for the k = 9 dipole mode with a period of 3 560 s 
of model L. In the homogeneous envelope model this mode 
is overstable but it is damped in the inhomogeneous one. 
Figure 11 displays in the lower panel the cumulative work; 
the solution for the homogeneous case is delineated with 
the full and the inhomogeneous one with the dashed line. 
Both work integral curves arc independently normalized to 
unity at the surface. Obviously, the hydrogen partial ion- 
ization zone contributes to driving in the homogeneous but 
damps in the inhomogeneous envelope, and it is the action 
of this region which decides over stability or overst ability. 
We attribute the swapping of the work contribution of the 
H-ionization zone to the influence of the different envelope 
structures. Notice that the kt runs are very close in both 
cases. Therefore, it is the phase relation between density and 
pressure perturbation at the H-ionization zone which cause 
the difference. Notice once more that that He drainage ex- 
tends only to log P w 5. 

3.4 Dependence on spherical degree 

Observed frequency spectra of roAp stars hint at the possi- 
bility that not only dipole modes are involved in their light 
variability. In some cases - such as HR 1217, HD 119027, 
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Figure 10. The run of Sj as a function of period for model L with 
weak temperature inversion. The full line shows the results for the 
inhomogeneous envelope. The dashed line depicts the behavior of 
oscillation modes in the model with homogeneous outer layers. In 
the inhomogeneous case, non-successive overstable radial orders 
exist. 
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Figure 11. Effect of He drainage in the outermost part of the 
envelope on the work integral of £ = 1, A; = 9 mode. The full line 
depicts the solution on the homogeneous envelope. The dashed 
line stands for the one on the chemically inhomogeneous envelope. 

HD 203932 - the frequency spacings are not even enough to 

be explainable with I — 1 modes only (Kurtz 1995). Also 
from the theoretical side, assuming now that a classical k 
mechanism drives the modes, we expect not only £ = 1 
modes to be excited. 

Figure 12 shows results for the inhomogeneous model L 
with weak temperature inversion. We plotted directly the 
imaginary parts of the eigenfrequencies for radial orders 
33 to 37. The highest order which becomes overstable at 
£ = 1 is 36. Lower orders tend to have their maximum 
driving at higher degrees. From the asymptotic relation 
i^ke oc j^o(fc + £/2 + e) we deduce that the maximum driving 
occurs always at similar period, irrespective of the spheri- 
cal degree. Our computations showed that a mix of different 
spherical degrees should contribute to the roAp variability. 
As it is monitored mainly photometrically at the moment, 
spherical degrees higher than 3 are not very relevant. Adding 
up all the \ow-£ (0 to 3) modes and the different radial or- 
ders, we found that usually less than 12 modes (in terms of 



Figure 12. Imaginary parts of high-order p modes of model L as 
a function of spherical degree. 

the principal quantum numbers k and £) are simultaneously 
excited. 

3.5 Node distribution in the atmospheres 

For recent spectroscopic observations of a Cir, phase 
shifts in the temporal behavior of absorption-line features 
which probe different depths in the atmosphere suggested 
that a node is lying in the optically thin superficial lay- 
ers (Baldry et al. 1998). We analyzed our eigensolutions in 
this respect. Figure 13 shows the overstable dipole modes of 
model H. The upper panel depicts displacement oigcnfunc- 
tions. They do not shift significantly for successive radial 
orders. In all three CEises displayed, a node lies at very low 
optical depths and a second one close to r = 2/3. The situa^ 
tion is different for the temperature perturbation for which 
the real part of the eigcnfunction is shown in the lower panel 
of Fig. 13. The A; = 35 mode shows a node at small opti- 
cal depths, whereas the next two lower-order modes have no 
node at r < 2/3. The reason for the temperature pertur- 
bation to have a stronger variation than the displacement 
function is the very superficial H/He ionization zone. This 
region causes, due to a strong entropy perturbation, an as- 
sociated strong temperature variation. 

We emphasize that our treatment of the outer regions, 
in particular of the optically thin regions are far from realis- 
tic. Therefore, the results in Fig. 13 should at best be taken 
as indications that indeed nodes can be expected in the at- 
mospheres of roAp stars. Quantitative aspects will certainly 
change when studies concerning this issue are performed 
with the correct perturbed transport equations. Only then 
can we expect to exploit this property of roAp stars for at- 
mospheric diagnosis. 

4 DISCUSSION 

Nonadiabatic stability analyses applied to evolutionary 
models with homogeneous Pop I composition revealed that 
the K mechanism of H/He ionization can destabilize short- 
period roAp-type modes. We found the short-period insta- 
bilities to be restricted to the 1.5 and 1.6 M© sequences. 
The rough location of the blue edge for the high-order roAp 
oscillations lies, according to our preliminary computations, 
at about log Teff = 3.86. Both, main-sequence and eaxly sub- 
giant models can become roAp pulsators. In contrast to a 
recent mapping of roAp stars onto the HR plane (North et 
al. 1997), we find our overstable models to have somewhat 
lower masses and to have lower effective temperatures than 
what the North et al.(1997) data imply. On the other hand, 
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Figure 13. Spatial behavior of eigcnfunction in the outermost 
layers for the 3 overstable high-order dipole modes of model H. 



wc expect the blue edge to shift to higher temperatures if 
the more helium (than what we assumed in our computa- 
tions) is depleted from the most superficial stellar layers. In 
this case the H driving in the H/Hel partial ionization zone 
should get more efficient. This has to be proven, however, 
with explicit computations in the future. 

In contrast to 5 Set stars, H/He ionization dominates 
the driving of the roAp modes. Since the H/He ionization 
zone lies at lower temperatures than the He"*" one, we expect 
the blue edge of the roAp instability domain to be somewhat 
cooler than that of 6 Set stars (cf. Fig. 1). 

In all cases, the number of overstable oscillation modes 
we counted in our computations was small. Not more than 
5 overstable modes of a given spherical degree were encoun- 
tered. This is still larger than what is observed in some roAp 
stars so that mode selection remains a problem to be ad- 
dressed. In contrast to 5 Set stars, however, the number of 
overstable modes is considerably smaller in roAp models. 
From Fig. 5 we see that at the base of the kt gradient in- 
duced by H/He partial ionization the thermal scale levels off 
at about 1 000 s. Only periods with shorter periods can be 
potentially excited in this region. As the overtones are al- 
ready high, the radiative dissipation is also high. Therefore, 
it seems that only for a few modes with high radial order can 
H/He driving overcome the radiative dissipation anymore. 

A necessary prerequisite for short-period oscillations to 
be excited was the hypothesis of a temperature inversion 
in the inner atmosphere of the models. The temperature 
inversion which raises the superficial critical frequency suf- 
ficiently to reflect roAp-type modes need not be very high. 
Depending on the steepness, one to three thousand K are 
enough. Another effect of the steepness of the inversion is 
that it influences the number of excited modes. 

The evolutionary stage of roAp stars can be anywhere 
between the main sequence and the beginning of the sub- 
giant stage - observed frequencies do allow for this band- 
width of dating. In the most advanced phases of evolution 
the core is no longer convective. Therefore, no magnetic dy- 
namo could be operative anymore; most of the star - except 
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for the narrow ionization regions - is radiative. If the cychc 
secular frequency variations of roAp stars were interpreted 
to be due to some magnetic activity cycle then the ohmic 
decay time of the stellar magnetic field would limit the max- 
imum ago of such roAp stars. At present, the origin of the 
frequency variations as observed by now in HD 134214 and 
HR 3831 (Kurtz et al. 1994, Kreidl et al. 1994, Kurtz et al. 
1997) is far from clear. 

EYom our computations we cannot unambiguously at- 
tribute an evolutionary stage to the range of periods of over- 
stable oscillation modes. From Fig. 7 we see a tendency of 
longer periods to be excited at later evolutionary stages in 
the I.SMq models. This finding could, however, not be re- 
produced in the 1.6Mq models. 

Introducing an atmospheric temperature inversion in 
the models eventually produced the long elusive overstable 
roAp-type modes. However, in the same models we found 
roAp and 6 Set modes excited - this contradicted observa- 
tional evidence. Therefore, we argued for a He drainage in 
these superficial regions where magnetic pressure dominates 
over gas pressure. When assuming a kilo-Gauss dipole field 
we found a mass fraction of only a few times 10~^M« to be 
magnetically dominated and hence drained from He. In such 
models we got rid of the S Set modes but could retain the 
roAp ones. As only the the H/He ionization was deficient 
of He, it is evident that H ionization dominates the excita- 
tion of roAp modes (see also the work integrals in Fig. 8). 
Notice that the He^ partial ionization is not modified in 
these models, but still we got rid of the S Set-type modes. 
This contradicts previous reasoning (Cox et al. 1979). Our 
computations of inhomogeneous models showed that for 5 
Set- type modes the H/Hc ionization zone has a damping 
influence which even over-compensates the driving action 
of the He+ region. Furthermore, observations support our 
picture which emphasizes the role of the magnetic field for 
the presence of high-frequency (roAp-like) oscillation modes 
which are not found in 6 Set stars such as in FG Vir (Breger 
et al. 1996). 

In some of the inhomogeneous models we found excited 
non-successive overtones. Such a behavior is reminiscent of 
the frequency spectrum as found for example in HD 217522 
(Kurtz 1995). 

A major problem with whole stability analysis as pre- 
sented here remains: The driving region of the roAp modes 
reaches into the region where the magnetic pressure dom- 
inates the star, therefore we must expect a non-vanishing 
influence of the perturbed magnetic field. In other words, 
magneto-acoustic modes might be excited which carry away 
energy (Roberts & Soward 1983, Dziembowski & Goode 
1996). If such an interaction is efficient enough it can damp 
out the pure acoustic modes of which we were talking here. 
However, at the moment we are not in a position to do such a 
nonadiabatic analysis quantitatively. Dziembowski & Goode 
(1996) presented estimates of such an interaction in the adi- 
abatic limit. From their results we expect that only the very 
strongest excited modes might survive in the correct hydro- 
magnetic treatment. 

Martinez (1996) discussed an unusual, possible roAp 
star - HD 75425 - which oscillates with 30-min period. This 
period is significantly longer than what is found for roAp 
stars and somewhat short for S Set stars. Referring to Fig. 7 
we find that such long periods could indeed be excited. For 



them the action of the He"*" is important, however. If we 
interpret HD 75425 in our framework, it should either be 
rather young so that sedimentation was not yet efficient 
enough or the stabilisation of the flow by the magnetic field 
is not sufficiently strong. The latter alternative would re- 
quire a magnetic field which is weaker than what is usually 
found in other roAp stars. It is, however, not yet clear why 
no high-order modes are present in HD 75425. Either they 
have not yet been seen and they are actually hidden (at 
even lower amplitude) in the long-period data or the hy- 
pothesized temperature inversion is not strong enough to 
trap the high-order p modes. 

It appears highly rewarding to search for weakly mag- 
netic A-type stars to see if it is possible at all to have both 
period domains (cf. Figs. 6 and 7) excited in any such stars. 
To extend our speculations further in this direction we men- 
tion Am - Fm stars (Kurtz et al. 1995, Kurtz 1998). These 
mildly peculiar cool main-sequence stars seem to be always 
binaries. Sedimentation of elements due to suppressed rota- 
tional velocities might eventually suppress 5 Set pulsations. 
Since the chemical peculiarity is not associated with a mag- 
netic field these stars cannot build up a sufficiently high 
critical frequency to support roAp-type modes. Therefore, 
we do not expect Am - Fm stars to be seen as short-period 
pulsators; on the other hand, if He drainage is not strong 
yet (i.e. if they are young enough) they might well be 5 Set 
type pulsators (cf. Cox et al. 1979). 

As mentioned before, we were able to enhance the crit- 
ical frequency at the stellar surface by ad hoc postulating 
a temperature inversion in the optically thin regions of the 
stellar atmosphere. Physically, this is equivalent to assum- 
ing a chromosphere in these stars. The spectral types of 
the roAp stars are just at the border where chromospheric 
evidence starts to be observed. According to prevailing the- 
oretical ideas, the roAp stars provide already the necessary 
ingredients for chromospheric heating: a superficial convec- 
tion zone and strong magnetic fields which might provide 
the necessary acoustic flux and possible magneto-acoustic 
waves, respectively. Shore et al. (1987) found, however, no 
signs of chromospheric emission in roAp stars. For further 
progress it would be highly rewarding to find observational 
evidence for or against the temperature-inversion hypothe- 
sis. 

The appendix discussed the approach of nonadiabatic 
frequencies to what is expected from adiabatic asymptotics. 
We found that caution is required when using high-frequency 

oscillation data to infer stellar parameters. Even if the ex- 
citation rates of high-order modes are low, the deviations 
of frequency spacings from the adiabatically predicted ones 
can be up to about 8% (see Fig. 15). Furthermore, this be- 
havior is expressedly non-monotonic so that an easy correc- 
tion cannot be applied. The magnitude of nonadiabatic fre- 
quency shifts in roAp stars is not negligible compared with 
the magnetic-field induced frequency perturbations (Dziem- 
bowski & Goode 1996) . Therefore, asteroseismology of roAp 
stars is quite a challenge for pulsation theory. To finish, we 
feel that a warning of simple application of adiabatic asymp- 
totics to high-order modes is also in place for solar seismol- 
ogy. 
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Figure 14. Frequency behavior for model B. The top diagram 
shows the imaginary parts of the nonadiabatic eigenvalues. The 
values for degrees 1 < ^ < 4 are plotted on top of each other. The 
middle panel shows the fractional deviation of the nonadiabatic 
frequency spacing from the asymptotic value. The lowest panel 
shows the same deviation computed from the adiabatic frequen- 
cies. The line coding is as following: £ = 1: full line; £ = 2: dotted 
line; £ = 3: dashed line; £ = 4: dash-dotted line. 
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APPENDIX 

We discuss concisely the effects of nonadiabaticity on 
the oscillation-frequency spectrum and compare it with adi- 
abatic results. This is mainly of interest for those aiming at 
using roAp frequency data for seismological deductions. 

For high radial-order modes of low degree, the frequency 
separation of equal-degree modes is usually rather close to 
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Figure 15. Relative deviation of the frequency spacing from the 
asymptotic value for model K. 



with the adiabatic sound speed given by Cad = y'FiP/p (cf. 
Sect. 3.4). The value of uq corresponds to the frequency of 
an acoustic wave to travel across a star. 

In the following we investigate how quickly the fre- 
quency separation approaches the asymptotic limit which is 
frequently referred to for asteroseismic purposes. To quantify 
the relative deviation of the actually computed frequency 
spacing from the expected one (z/o), we introduce 



^ = 1. 

1^0 



I'k + l 



(A2) 



1=2 

VO Jo Cad 



(Al) 



Figure 14 shows extensive results for model B for which 
we computed oscillation modes with spherical degree 1 (full 
hue), 2 (dotted line), 3 (dashed line), and 4 (dash-dotted 
line) for many radial orders. Computations were performed 
adiabatically as well as non-adiabatically. The bottom panel 
of Fig. 14 displays the relative deviation dv/vo of the adia- 
batic frequency spacing as a function of the oscillation fre- 
quency. Above about 1.3 mHz the deviation from the asymp- 
totic value (51.46 /^Hz) drops below 2 % for all degrees. 
In the range between 2.25 and 2.5 mHz, where many roAp 
modes are observed, the adiabatic mode separation agrees 
with the asymptotic value on the level of 0.1 to 1 %, de- 
pending on the degree. This frequency range corresponds to 
radial orders exceeding 40. 

To sec if the convergence properties of the nonadiabatic 
eigenfrequencies are comparable with the adiabatic ones we 
plot the relative deviation from the asymptotic value (adia- 
batically computed) in the middle panel of Fig. 14. The most 
obvious discrepancy is seen around 1 and 1.3 mHz. There 
we find a relative deviation of the spacing of 40 % and 10 %, 
respectively. These two features arc clearly associated with 
the local maxima in the damping rates ai which are shown 
in the top panel of Fig. 14. In the frequency range between 
1.5 to 2 mHz, the frequency separation relaxes to within a 
few percent of the asymptotic one. For higher frequencies, 
however, they start to diverge again in accordance with the 
rise of the damping rates. 
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What relevance does this behavior have for roAp 
modes? After all, the modes which we used for the above 
demonstration were all damped. We consider now the rel- 
ative deviation Av/pq of the nonadiabatic oscillation fre- 
quency separation from the the adiabatic asymptotic value 
iov i — 1 of model K. This model has 4 overstable modes 
around 500 s period. The imaginary parts are rather small 
but still the frequency spacing reacts sensitively to the cri 
variation. Figure 15 shows that at 1.5 mHz - at the position 
of a marginally overstable oscillation mode for which the 
damping rate is close to zero - the frequency spacing devi- 
ates by about 8 % from the adiabatically expected one. lu 
the overstable period domain, the frequency spacings change 
from 3 % longer to 8 % shorter than the predicted asymp- 
totic frequency spacing. For seismological inferences such 
differences can be relevant. 

To finish, we emphasize that, if high accuracy is crucial, 
nonadiabaticity effects are not negligible in deriving stellar 
parameters from high-frequency oscillation data. This con- 
clusion might not only apply to roAp stars but also to the 
high-frequency domain of solar oscillations. 
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